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Abstract 

Background: This study aims to research two areas, one with a resistant and the other with a susceptible profile 
of An. gambiae to deltamethrin in the region of Plateau (southern Benin). In each area, eight localities were sought. 
Both areas were needed for the assessment of the impact of malaria vector resistance to pyrethroids on the 
effectiveness of Long Lasting Insecticidal Nets (LLINs). The susceptible area of An. gambiae to deltamethrin was 
used as a control. 

Methods: In total, 1 19 localities in the region of Plateau were screened by sampling An. gambiae s.l larvae. Female 
mosquitoes resulting from these larvae were exposed to 0.05% deltamethrin following WHO standards. PCR was 
used to identify species and molecular forms of the dead and alive mosquitoes. Finally, we identified kdr mutations 
(1014 F and1014S) using the HOLA technique. 

Results: Fifty-six out of 1 19 prospected localities tested positive for Anopheles gambae s.l breeding sites. The results 
showed that An. gambiae was resistant to deltamethrin in 39 localities and susceptible in only 2 localities; resistance 
to deltamethrin was suspected in 15 localities. The HOLA technique confirmed the presence of kdr 1014 F mutation 
and the absence of kdr 1014S mutation. The kdr 1014 F mutation was found in both M and S molecular forms at 
relatively high frequencies therefore confirming the susceptibility tests. 

Conclusion: We were unable to identify the eight susceptible areas due to the overall resistance of An. gambiae to 
deltamethrin in the region of Plateau. To implement the study, we kept two areas, one with high resistance (R +++ ) 
and the other with low resistance (R + ) of An. gambiae to deltamethrin. 
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Background 

Benin, like many other African countries, has based its 
vector control strategy on two major interventions: uni- 
versal access to Long Lasting Insecticidal Nets (LLINs) 
and Indoor Residual Spray (IRS). Thus, in late July 2011, 
a large scale distribution of LLINs had been conducted by 
Benin National Malaria Control Program (NMCP) with 
the support of the World Bank, the US President's Malaria 
Initiative (PMI) and the World Health Organization 
(WHO). Another distribution campaign had already been 
conducted in October 2007 to ensure a partial and select- 
ive coverage in two types of LLINs, PermaNet 8 and 
OlysetNet 8 respectively impregnated with deltamethrin 
and permethrin. A national survey assessing the 2007 
campaign revealed that 56.3% of children under 5 years 
and 54.8% of pregnant women slept under a LLIN the 
night before the survey [1]. The current study carried out 
between April and early July 2011 had taken place be- 
tween the two aforementioned distribution campaigns. 

Insecticide resistance has become widely distributed in 
Western [2-5], Eastern [6], Central [7] and Southern 
Africa [8]. This could be a serious obstacle to the effect- 
iveness of LLINs. This raises an important question 
within National Malaria Control Programmes (NMCPs): 
should we continue to promote LLINs? The question 
has been explored by N'Guessan et al. [9] who demon- 
strated a decrease in the effectiveness of LLINs and 
lambdacyhalothrin IRS in an area of high resistance of 
Anopheles gambiae in Southern Benin. They compared 
the effectiveness of mosquito nets in two areas, one area 
with resistant An. gambiae populations to pyrethroids and 
one other where An. gambiae is susceptible to pyrethroids. 
However, the study was conducted on an experimental 
scale using experimental huts. It is therefore difficult to 
predict what would happen at a community level. 

For that reason, a large scale study investigating the 
impact of malaria vector resistance to pyrethroids on the 
effectiveness of LLINs was initiated in late July 2011. In 
the context of universal access to LLINs, the implemen- 
tation of such an impact study requires the identification 
of two areas: one area where An. gambiae was highly 
resistant to pyrethroids, and one other where An. gambiae 
was still susceptible to pyrethroids. The latter area would 
serve as a control area. The main objective of this article 
is the determination of these two areas. This will allow 
the comparison of the effectiveness of LLINs between 
the two areas. 

The results of a previous study conducted in forty 
districts in Southern Benin revealed that the entire re- 
gion is covered by pyrethroids resistant populations of 
An. gambiae [5]. However, a few pockets of susceptible 
populations of An. gambiae has been reported in two 
districts (Ifangni and Adja-were) located in the region of 
Plateau [5]. This is why we selected this region for the 



implementation of our impact study. However, we de- 
cided to randomly implement the study in 16 localities 
(8 localities constituting the resistant area and 8 
others representing the susceptible area). We then 
multiplied larval prospections and susceptibility tests 
in many localities across four districts within the re- 
gion of Plateau. Unfortunately, we were unable to identify 
the 8 localities of high susceptibility of An. gambiae as 
previously expected. We then decided to replace the 
susceptible area by an area where An. gambiae has a low 
resistance [low resistant area (R + area) vs high resistant 
(R +++ area)]. 

Regarding the susceptibility tests, we had the choice 
between deltamethrin and permethrin, two insecticides 
used to impregnate PermaNet 8 and OlysetNet 8 . We ended 
up choosing deltamethrin. This choice has been motivated 
by evidence suggesting a wide distribution of high vector 
resistance to permethrin in Southern Benin [5]. Under 
these conditions, we believed that the probability to iden- 
tify localities where An. gambiae is susceptible to per- 
methrin is low. 

From April to early July 2011, we conducted several 
field missions across four districts in the region of Plateau. 
Results from the susceptibility tests helped classify the 
localities between the two resistant (R + and R +++ ) areas. In 
both R + and R +++ areas, a study of malaria transmission 
and vector behavior had been implemented. 

Methods 

Study area 

The study was conducted in the region of Plateau in 
the south-east of Benin and specifically in the dis- 
tricts of Ifangni, Sakete, Pobe and Ketou. This region 
has an area of 3,264 km , with a total population of 
407,116 inhabitants [10]. Fifty-six localities which tested 
positive for An. gambiae larvae are mostly rural and 
represent 31.3% (56/179) of the localities of the four 
districts to which they belong, which total area is 
2,849 km 2 [10]. The climate in the Plateau region is 
Guinean with two rainy seasons and two dry seasons. 
The region records an annual rainfall between 
800 mm and 1200 mm in its western part and be- 
tween 1000 mm and 1400 mm in its eastern part. 
The cropping system is characterized by the practice 
of two annual growing seasons associated with rainfall 
patterns. Corn crop is predominant. The whole area 
is strewn with swamps. These swamps are used for 
the production of off-season crops and the installation of 
ground farms for various species. Different economic 
activities take place in the region of Plateau because of the 
opportunities offered by the natural environment but also 
its closeness to Nigeria, which has a polarizing action in 
the region. The region of Plateau abounds with farmers 
and traders. 
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Collection and breeding of larvae and pupae of 
Anopheles mosquitoes 

From April to early July 2011, 119 localities were pros- 
pected in all four districts, of which 56 tested positive 
for Anopheles breeding sites. Localities which tested 
positive were geo-referenced using a Global Positionning 
System (GPS). 

Larvae and pupae of An. gambiae were collected in collec- 
tions of water using the "dipping" technique [11]. This tech- 
nique consists of collecting larvae and pupae, using a simple 
ladle, from positive breeding sites. Larvae and pupae were 
kept separately in labeled bins and taken to the insectary 
of the "Centre de Recherche Entomologique de Cotonou" 
(CREC) for rearing. After emergence, adults of 2 to 5 days 
old were used for susceptibility testing in the laboratory. 

Susceptibility tests 

Susceptibility tests were performed according to the 
WHO susceptibility tube-test with unfed female Anoph- 
eles, aged 2 to 5 days. These tests were performed with 
paper impregnated with deltamethrin at the diagnostic 
dose of 0.05%. 

Batches of 25 female mosquitoes aged 2 to 5 days were 
added to each tube carpeted with deltamethrin impreg- 
nated paper for 60 minutes. Batches exposed to untreated 
papers were used as control. Mosquitoes of the susceptible 
Kisumu strain were exposed to deltamethrin (0.05%) trea- 
ted filter papers. The number of field mosquitoes knocked 
down as a result of the insecticide effect was recor- 
ded every ten minutes during the exposure time period. 
After 60 minutes of exposure, mosquitoes were trans- 
ferred to the observation tubes and fed with a 10% honey 
solution and kept under observation for 24 hours. At the 
end of the observation time, mortality rates were deter- 
mined. These rates were construed in accordance with the 
recommended criteria by WHO [12]. The resistance status 
was determined based on the following criteria: 

- Mortality > 97%: susceptible Anopheles population. 

- Mortality 80 - 97%: suspected resistance in the 
Anopheles population. 

- Mortality < 80%: resistant Anopheles population. 

Given the small number of localities where An. gambiae 
was susceptible to deltamethrin, we decided to define R + 
localities as all localities where Anopheles mortality was 
greater than or equal to 80%. All localities with mortal- 
ity rates below 80% were considered highly resistant 
(R +++ localities). 

We did not choose the 80% cutoff at random to separate 
both areas. This criterion was chosen according to the 
WHO standard to identify highly resistant mosquitoes 
[12]. We considered suspected resistant mosquitoes and 
susceptible ones [12] as lowly resistant. 



After susceptibility testing, mosquitoes were kept on 
silicagel at -20°C for molecular characterization. 

Selection of the localities into R + and R +++ areas 
After testing for susceptibility, we randomly selected the 
sixteen localities for the implementation of our impact 
study at both defined resistance levels (R + and R +++ ) in 
order to avoid selection bias. 

Molecular characterization of An. gambiae, PCR species, 
molecular forms and kdr 1014 F and 1014S 

For each locality, 22-49 females Anopheles were analyzed 
by PCR according to the protocol described by Scott et al. 
[13]. Identification of Anopheles species was made accord- 
ing to the protocol described by Favia et al. [14]. For each 
locality, DNA extracts from An. gambiae allowed us to 
look at the simultaneous presence of kdr mutation 1014 F 
and 1014S as recommended by Lynd et al. [15]. 

Mapping of the resistance 

The geographical coordinates of the localities where larvae 
were collected from were recorded by GPS and projected 
onto a map of the region of Plateau. The three levels of 
susceptibility of Anopheles to deltamethrin (0.05%) accor- 
ding to the criteria of the WHO [12] (sensitive, resistant 
and suspected resistance) were shown on the map. Over- 
all, 56 localities were mapped. 

Statistical analysis 

The calculation of mortality rates was performed using 
MS Excel spreadsheet. The chi-square test of compari- 
son of proportions was used to compare mortality rates 
within localities of the same district and within the two 
defined areas (R + and R +++ ). Knockdown times (KdT50 
and KdT95) were determined by logistic regression with 
probit link. We used the delta-method to determine the 
95% Confidence intervals of the KdTs for the two areas. 
The variation in kdr frequency between the two areas, 
and between both M and S forms of An. gambiae was 
assessed via a logistic regression [16]. The same statis- 
tical method was used to assess the spatial variation in 
resistance levels within districts. The analysis of deviance 
penalized by the dispersion parameter [17] was used to 
assess the relevance of the variability. 

All analyses were performed with R-2.15.2 statistical 
software [18]. 

Results 

Deltamethrin-induced mortalities in 56 localities 

The An. gambiae Kisumu reference strain was suscep- 
tible to deltamethrin (0.05%), showing 100% mortality. 

Mortality rates obtained after exposing An. gambiae 
wild populations to deltamethrin (0.05%) varied between 
20% and 100%. 
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In Ifangni, mortality rates ranged from 20% in Gblo-gblo 
and Djegou-djegui to 100% in Zihan. Out of 19 localities, 
An. gambiae was found susceptible only in 2 localities 
(Zihan and Ko-Ai'djedo). Suspected resistance was observed 
in 10 localities (Table 1). 

In Sakete, An. gambiae was not susceptible. In Iwa'i, 
Idi and Djohounkolle suspected resistance was noted 
(Table 1). 

In the district of Pobe, mortality of An. gambiae mosqui- 
toes ranged from 45% in Obanigbe to 85% in Agbarou. An. 

Table 1 Resistance status of An. gambiae s.l populations 
to deltamethrin in the localities of Ifangni and Sakete 



Districts Localities 


Total 

(N) 


KdT50 
(min) 


KdT95 
(min) 


Dead 
(N) 


Mortality 

(%) 


Ifongni Lokossa* 


44 


38.9 


74.0 


36 


82 a 


Ko-dogba* 


48 


36.1 


68,6 


42 


87 ad 


Zihan* 


50 


60.0 


114.0 


50 


100 


Baoudjo 


89 


40.5 


77.0 


72 


81 a 


Gbedji 


88 


100.0 


190.0 


68 


77 a 


Igolo 


22 


111.1 


211.1 


1/ 


77 a 


Araromi* 


39 


41.6 


79.2 


33 


85 ad 


Itassoumba* 


66 


39.2 


74.5 


59 


89 ad 


Ko-koumolou* 


48 


32.6 


62.0 


40 


83 a 


Itakpako* 


131 


36.6 


69.5 


105 


80 a 


Banigbe Centre 


97 


42.3 


80.3 


32 


33 bc 


Tchaada 


22 


93.8 


178.1 


/ 


32 bc 


Daagbe* 


40 


31.6 


60.0 


33 


82 a 


Akadja 


1 23 


38.5 


73.1 


62 


JU 


Gblo-gblo 


45 


103.4 


196.5 


9 


20 c 


Ko-aTdjedo* 


92 


31.2 


59.4 


90 


98 d 


Zoungodo* 


62 


43.5 


82.6 


52 


84 a 


Djegoun-djegui 


40 


75.0 


142.5 


8 


20 c 


Ketougbekon* 


147 


35.7 


67.8 


121 


82 b 


Sakete llakofadji 


44 


31.6 


60.0 


30 


68 a 


Dagbao 


45 


32.2 


61.3 


33 


73 a 


Itadjebou 


34 


42.8 


81.4 


21 


62 abc 


Djohounkolle* 


73 


39.4 


75.0 


59 


81 a 


Igbo-abikou 


65 


115.4 


219.2 


25 


38 b 


Igbola 


38 


96.7 


183.8 


19 


50 ab 


Alabansa 


30 


150.0 


285.0 


12 


40 bc 


IwaT* 


43 


37.0 


70,3 


37 


86 a 


Ikemon 


86 


71.4 


135.7 


40 


46 bc 


Idiagbola 


40 


52.6 


100.0 


1/ 


43 bc 


Yoko Centre 


40 


31 


58.7 


30 


75 a 


Idi* 


31 


62.5 


118.7 


30 


97 d 



gambiae proved resistant in all localities except Agbarou 
where suspected resistance was observed (Table 2). 

In Ketou, the lowest mortality rate of An. gambiae (67%) 
was obtained in Alakouta versus 83% in Adjozounme. Out 
of a total of 15 localities evaluated, suspected resistance 
was noted only in Adjozounme, while full resistance was 
observed in all other localities (Table 2). 

The lowest KdT50 (31 minutes) was obtained in Yoko 
Centre (district of Sakete) versus 375 minutes in Illekpa 
(district of Pobe) (Tables 1 and 2). 

Overall, across the 56 localities of the four districts, An. 
gambiae was resistant to deltamethrin in 39 localities. 
However, susceptibility was noted in An. gambiae mosqui- 
toes from 2 localities of the district of Ifangni. Resistance 
of An. gambiae to deltamethrin (0.05%) was suspected in 
15 localities (Tables 1 and 2). 



Table 2 Resistance status of An. gambiae s.l populations 
to deltamethrin in the localities of Pobe and Ketou 



Districts 


Localities 


Total 

(N) 


KdT50 
(min) 


KdT95 
(min) 


Dead 
(N) 


Mortality 


Pobe 


Igbo-okpa 


79 


85.7 


162.8 


58 


73 3C 




Okeita 


1 01 


47.6 


90.5 


/8 


77 ac 




Obanigbe 


40 


57.7 


109.6 


18 


45 b 




Agbarou* 


55 


66.6 


126.6 


4/ 


85 a 




Issaba 


111 


100 


190 


68 


61 c 




Okoofi 2 


54 


96.7 


183.8 


24 


44 b 




Illekpa 


39 


375 


712.5 


26 


67 c 




Osoumou 2 


48 


53.6 


101.8 


34 


71 c 




Osoumou 1 


70 


142.9 


271.4 


44 


63 c 




Onigbolo 


90 


52.6 


100 


66 


73 ac 


Ketou 


Alakouta 


36 


68.2 


129.5 


24 


67 abc 




Okpometa 


5/ 


81.1 


154 


23 


41 a 




Mowodani 


52 


68.2 


129.5 


31 


60 ab 




Idena 2 


179 


39 


74 


81 


45 a 




Kpankoun 


70 


65.2 


124 


37 


53 ab 




Okeola 


49 


66.6 


126.6 


24 


48 a 




Odokoto 


95 


33.7 


64 


51 


54 ab 




Igui-olou 


100 


49.2 


93.4 


72 


72 bc 




Adjozoume* 


100 


33.3 


63.3 


83 


83 c 




Idena 3 


80 


68.6 


129.5 


39 


49= 




Oloumou 


94 


50.8 


96.6 


48 


51 ab 




Kouhoudou 


60 


47.6 


90.5 


37 


62 ab 




Obatedo 


95 


43.5 


82.6 


74 


78 bc 




Atchoubi 1 


100 


33 


62.6 


// 


77 bc 




Atchoubi 2 


55 


61.2 


116.3 


35 


64 bc 



a, b, c, ab, ac, morta |j t y rates w j tn tne different superscript in the same district 
are statistically different (p < 0.05). 

*The localities with an asterisk are lowly resistant and those without an asterisk 
are highly resistant. 



a, b, c, ab, ac-p^ mor taMty rates with the different superscript in the same district 
are statistically different (p < 0.05). 

The localities with an asterisk are lowly resistant and those without an asterisk 
are highly resistant. 
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Categorization of localities as R + and R +++ across the four 
districts 

Overall, out of the 56 localities, 16 were categorized as R + 
and 40 as R +++ based on the mortality rates observed 
(Tables 1 and 2). Out of the 16 R + localities, 11 including 
the two susceptible ones were located in the district of 
Ifangni, 3 in Sakete, 1 in Pobe and 1 in Ketou (Tables 1 
and 2). Localities classified as R +++ were distributed as 
follows: 8 in Ifangni, 9 in Sakete, 9 in Pobe and 14 in Ketou 
(Tables 1 and 2). 

Based on these results, 8 localities (Itakpako, Araromi, 
Ko-koumolou, Djohounkolle, Ko-Aidjedo, Lokossa, 
Ketougbekon and Iwa'i) were randomly selected in the R + 
area and 8 other (Mowodani, Banigbe, Okoofi 2, Ikemon, 
Akadja, Idena 2, Igbola and Tchaada) in the R +++ area for 
the implementation of our impact assessment study. 

Mortality rates were 84% (516/617) in the R + area 
made of 8 R + localities versus 45% (296/651) for the R +++ 
area composed of 8 R +++ localities (p <0.0001). Statistical 
analysis also showed that the risk that mosquitoes die in 
contact with deltamethrin was six times higher in the R + 
area than in the R +++ area (OR = 0.15, 95% CI: 0.10-0.23, 
p < 0.0001). In addition, the KdT95 was 67 minutes (95% 
CI: 52.80-86.21) in R + area versus 93.4 minutes (95% CI: 
87.06-113.93) in the R +++ area (p < 0.05). 

Molecular characterization of mosquitoes from 16 selected 
localities (8 R + and 8 R +++ ) 

Identification of species and forms of An. gambiae complex 

After testing for susceptibility, Polymerase Chain Reaction 
(PCR) was used to identify species of An. gambiae s.l 



complex, and the identification of the different forms of 
each An. gambiae s.s species. This characterization was 
carried out on the collected mosquitoes from the 16 
selected localities. 

The results show that out of the 506 specimens of An. 
gambiae s.l evaluated in all 16 localities, only An. gam- 
biae s.s was found in all localities. Both molecular forms 
M and S of An. gambiae s.s were found in sympatry in 
all the 16 localities. Overall, 66.8% (338/506) of the mos- 
quitoes belonged to the form M and 32.6% (165/506) to 
the form S. The proportion related to one or the other 
form varies depending on the locality (Figure 1). An. 
gambiae s.s form M was predominant in most localities 
except Lokossa and Ko-Aidjedo in Ifangni, Igbola in 
Sakete and Idena 2 in Ketou. Both forms were found in 
similar proportions in Tchaada (12 An. gambiae s.s form 
M out of 22). 

Some hybrid individuals M/S (0.6%) were also found in 
Itakpako (district of Ifangni) and Igbola (district of Sakete). 

Detection of the kdr mutation in the molecular forms 
(M and S) of An. gambiae s.s 

The 1014S kdr mutation originally from East Africa was 
absent in 497 mosquitoes tested. However, the 1014 F 
kdr mutation was detected in the four districts of the 
region of Plateau. The 1014 F kdr mutation was found 
in both M and S molecular forms of An. gambiae s.s but 
at variable frequencies depending on the locality. Indeed, 
the frequency of the kdr mutation varies between 0.50 
and 0.93 in the form M and between 0.45 and 0.92 in 
the form S (Table 3). The small number of characterized 



Tchaada 
Okoofi 2 
Mowodani 
Ikemon 
Igbola I 
Idena 2 I 
Banigbe | 
Akadja 




13 

U Djohounkolle 
^ Araromi I 
Lokossa I 
Ko-koumolou 
Ko-Aidjedo 
Ketougbekon 
Iwai 
Itakpako 

20 40 60 80 

Frequency of molecular forms (%) 

Figure 1 Distribution of An. gambiae s.s. molecular forms per locality. 
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Table 3 Frequencies of 1014 F kdr mutation within 
An. gambiae s.s. molecular forms 



Localities/Areas 


Molecular 
forms 


N 


ss 


RS 


RR 


F 

[kdr) 


R + localities 














Itakpako 


M 


27 


0 


17 


10 


0.69 




S 


2 


0 


2 


0 


0.50 


Araromi 


M 


24 


3 


18 


3 


0.50 




S 


2 


0 


1 


1 


0.75 


Ko-koumolou 


M 


19 


1 


16 


2 


0.53 




S 


3 


0 


3 


0 


0.50 


Djohounkolle 


M 


26 


2 


10 


14 


0.73 




S 


1/ 


0 


3 


14 


0.91 


Ko-ATdjedo 


M 


5 


0 


4 


1 


0.60 




S 


28 


3 


16 


9 


0.61 


Lokossa 


M 


6 


0 


2 


4 


0.83 




S 


16 


0 


8 


8 


0.75 


Ketougbekon 


M 


15 


2 


6 


7 


0.67 




S 


4 


1 


3 


0 


0.38 


IwaT 


M 


30 


0 


4 


26 


0.93 




S 


3 


0 


1 


2 


0.83 


ft + area 


M + S 


227 


12 


114 


101 


0.70 


R +++ localities 














Mowodani 


M 


28 


0 


7 


21 


0.88 




S 


12 


1 


4 


7 


0.75 


Banigbe 


M 


28 


0 


10 


18 


0.82 




S 


5 


2 


1 


2 


0.50 


Okoofi2 


M 


43 


0 


12 


31 


0.86 




S 


6 


0 


1 


5 


0.92 


Ikemon 


M 


29 


0 


7 


22 


0.88 




S 


9 


0 


4 


5 


0.78 


Akadja 


M 


46 


1 


29 


16 


0.66 




S 


2 


0 


2 


0 


0.50 


Idena2 


M 


2 


0 


2 


0 


0.50 




S 


19 


2 


17 


0 


0.44 


Igbola 


M 


/ 


0 


6 


I 


0.57 




S 


12 


0 


8 


4 


0.66 


Tchaada 


M 


11 


0 


9 


2 


0.59 




S 


10 


0 


8 


2 


0.60 


fi +++ area 


M + S 


269 


6 


127 


136 


0.74 


Total [(R + and fi +++ J areas] 


M 


346 


9 


159 


178 


0.74 




S 


150 


9 


82 


59 


0.67 



mosquitoes in some localities did not allow for accurate 
estimation of the frequency of this mutation. These local- 
ities are Tchaada and Lokossa in the district of Ifangni and 
Igbola in the district of Sakete. 



By combining data from the two areas (R + and R +++ ), 
the kdr frequency was similar within both molecular 
forms: M [F{kdr) = 0.74] and S [F(kdr) = 0.67] with p = 0.17 
(Table 3). Furthermore, the aggregated data from the two 
molecular forms yielded a kdr gene frequency estimated at 
0.70 in the R + area versus 0.74 in the R +++ area. There is 
no difference between the frequency of the kdr mutation 
between the two areas (p = 0.13). 

Mapping of the resistance of An. gambiae to deltamethrin 

Various resistance levels to deltamethrin were observed 
based on the WHO criteria [12]. Figure 2 illustrates the 
mapping of the resistance of An. gambiae to deltameth- 
rin in the prospected localities. Overall, the map shows 
the profile of An. gambiae resistance to deltamethrin in 
the region of Plateau in 2011. A concentration of localities 
where An. gambiae was susceptible or suspected resist- 
ance to deltamethrin is distinguishable in the district of 
Ifangni. On the other hand, An. gambiae was resistant in 
most localities in the districts of Sakete, Pobe and Ketou. 

Table 4 shows the results of logistic regression analysis 
that was performed to assess differences in the level of 
resistance between the districts. In general, although the 
district of Ifangni was composed of localities where An. 
gambiae was susceptible or suspected resistant, it is clear 
from this analysis that the level of resistance has not 
been a significant variation within the four districts 
[p (LR-test) =0.22]. 

Discussion 

Our study shows a widespread resistance of An. gambiae 
to deltamethrin in the Plateau region, southeastern 
Benin. Taking into account the KdTs, the effect of delta- 
methrin on An. gambiae mosquitoes was not uniform 
across the various localities. Resistance with deltameth- 
rin was associated with relatively high kdr frequencies 
found in both M and S form of An. gambiae s.s. 

Out of the 56 localities explored, An. gambiae was 
found susceptible in only two localities (Ko-Aidjedo and 
Zihan). These results show a wide distribution of delta- 
methrin resistance as previously reported by Padonou 
et al. [19] in certain localities of Oueme, a department 
near Plateau region. The results also confirm the rapid 
expansion of pyrethroid resistance in natural populations 
of An. gambiae in Africa [20] and particularly in Benin 
[21-23]. However, the small numbers of mosquitoes ex- 
posed to deltamethrin in some localities did not allow 
accurate estimations of the mortality rates. This consti- 
tutes a major limitation for our study. 

Previous studies demonstrated that malaria vector re- 
sistance to pyrethroids might be related to an extensive 
and massive use of LLINs [24,25]. Insecticide molecules 
at the surface of LLINs might exert a lethal effect on 
susceptible mosquitoes, therefore selecting for resistant 
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Figure 2 Map showing the distribution of An. gambiae resistance to deltamethrin in the Plateau region. 
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Table 4 Spatial variation of the resistance level of An. gambiae to deltamethrin within districts 


Districts 


Total 


Dead (N) 


Dead (%) 


Coef 


OR 


95% CI 


p (Wald test) 


p (LR-test) 


Ifangni 


1293 


936 


72.39 


0.00 


1.00 






0.22 


Sakete 


569 


353 


62.04 


-0.56 


0.57 


[0.29-1.14] 


0.11 




Pobe 


687 


463 


67.39 


-0.24 


0.79 


[0.41-1.51] 


0.47 




Ketou 


1222 


736 


60.23 


-0.51 


0.60 


[0.35-1 .04] 


0.06 





mosquitoes that reproduce in natural populations. In the 
department of Oueme for instance, Padonou et al. [19] 
reported an increase in the resistance of An. gambiae to 
deltamethrin following a mass distribution campaign of 
LLINs. It is then possible that the resistance level has 
increased very quickly in the region of Plateau after the 
selective distribution campaign of LLINs in 2007. Several 
other factors such as the domestic use of pyrethroids [3] 
and water run off loaded with insecticide particles from 
the North of Benin where pyrethroids were massively 
and uncontrollably used to control cotton pests [26] could 
also explained the resistance of An. gambiae to deltameth- 
rin in the localities of the region of Plateau. The insecticides 
particles contained in the water could exert a selection 
pressure on the larvae of An. gambiae [27]. 

Our data surprisingly reveals a low resistance of An. 
gambiae to deltamethrin in the district of Ifangni. None- 
theless, we did not observe any spatial variation in the 
level of resistance across the four districts (p = 0.22). 
This suggests a similar selection pressure on An. gambiae 
s.l across the four districts. Therefore, other reasons, still 
unknown to our knowledge, may explain the low levels of 
resistance of Anopheles population to deltamethrin in the 
district of Ifangni. Further studies must be undertaken to 
understand the real causes of the low resistance of An. 
gambiae in this district. 

PCR analysis reveals that An. gambiae s.s was the only 
species of the An. gambiae s.l. complex encountered in 
all 16 selected localities. The absence of An. melas could 
be explained by the fact that the larvae were all collected 
from small freshwater pools. This is understandable 
since An. melas larvae are mainly found in brackish 
water ponds [28]. Similarly, An. arabiensis was not present 
even though it has already been reported in Central 
Benin [23]. 

Regarding the molecular forms, the relative dominance 
of one form over the other could be explained by the pres- 
ence of specific breeding sites to one or the other mole- 
cular forms [29]. Globally, in both areas (R + and R +++ ), 
66.8% of the mosquitoes were of M form, 32.6% of S form 
and 0.6% of the hybrid M/S. This contradicts Yadouleton 
et al. [5] who reported the absence of the form S of An. 
gambiae s.s in the region of Plateau. In fact the presence 
of the form S in significant proportions in 2011 in 
that region could be associated to the rainfall pattern 



and the rapid infiltration of water into the soil allowing 
the formation of favorable temporary breeding sites for 
the development of form S of An. gambiae s.s. 

Kdr mutation 1014 F has been found in M as well S 
form of An. gambiae s.s, but at variable frequencies 
depending on the locality. Djenontin et al. [30] reported 
a higher frequency of kdr mutation in the form S than in 
the form M of An. gambiae s.s in Ouidah-Kpomasse- 
Tori Bossito area between October and December 2007. 
The same trend has been reported by Diabate et al. [31] 
and Dabire et al. [32] in the Kou valley in Burkina-Faso 
and by Dabire et al. [33] in Guinea-Bissau. However, for 
our study, data from both areas (R + and R +++ ) revealed a 
similar kdr frequency between the two molecular forms 
(p = 0.17). This suggests that the selection pressure exer- 
ted on the two molecular forms in the natural environment 
did not differ significandy across the region. 

We did not find the 1014S kdr mutation originally 
from East Africa [34]. However, previous studies have in- 
dicated that this mutation was present in Benin [23] and 
was expanding to the north and center of Benin in An. 
gambiae. It is therefore important to extend the surveil- 
lance of Wrl014S in the region of Plateau. 

Although a significant difference in mortality rates was 
found between both areas (R + and R +++ ) (p < 0.0001), no 
difference was observed between the frequencies of kdr 
gene in both areas (p = 0.13). In addition, KdT95 was 
higher in R +++ area than in R + area (p < 0.05). These 
findings suggest the involvement of other resistance 
mechanisms in addition to the kdr mutation in the mos- 
quitoes from the R +++ area. 

Recent studies have also reported the involvement of 
certain metabolic enzymes in the resistance of malaria 
vectors to pyrethroids in Africa [35] and in several other 
regions in Benin [36,37]. For example, an overexpression 
of CYP6M2 and CYP6P3 genes, involved in the metabol- 
ism of pyrethroids, has been reported in resistant popu- 
lations of An. gambiae in Porto Novo [36]. Therefore, 
we believe that the role of metabolic resistance in the 
region of Plateau deserves further scrutiny. 

From the results recorded from the susceptibility tests, 
we were able to map the distribution of the resistance of 
An. gambiae to deltamethrin. Such a map was necessary 
since it provides a picture of the availability of localities 
where An. gambiae was still susceptible to deltamethrin. 
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In addition, this map could be used as an important tool 
to monitor the dynamics of the resistance of An. gambiae 
to pyrethroids. Moreover, given that our research activities 
were carried out in collaboration with Benin NMCP, this 
map should allow Benin NMCP to accordingly adapt its 
strategy of malaria prevention. 

Given the extensive use of LLINs inside houses at this 
time, it is possible that highly resistant mosquitoes would 
tend to feed inside more than low resistant Anopheles 
mosquitoes. The deterrence effect of LLINs is likely to be 
more effective on low resistant Anopheles population. If 
this was verified, the likelihood for highly resistant mos- 
quitoes to transmit malaria parasites would be higher than 
that of lowly resistant ones. If this assumption was true, it 
would confirm that vector resistance to insecticides is a 
major concern to the operational effectiveness of the 
LLINs distributed in 2011 [38]. 

Conclusion 

The results of this study show a wide distribution of the 
resistance of malaria vectors to deltamethrin in the re- 
gion of Plateau. An. gambiae was found susceptible in 
only two locations in the district Ifangni. Two popula- 
tions of An. gambiae s.s. were encountered in the region: 
An. gambiae s.s form M and An. gambiae s.s form S. Both 
molecular forms were resistant to pyrethroids. Besides the 
kdr mutation, our study suggested the involvement of 
other resistance mechanisms of An. gambiae to pyre- 
throids in the region of Plateau. 

Based on the resistance criteria we defined, more than 
70% of the 56 prospected localities were classified as R +++ 
localities. This confirms the rapid expansion of the resist- 
ance phenomenon across the region. Because of this ex- 
pansion, we were not able to find the eight susceptibility 
localities for the implementation of our impact study of 
the resistance of An. gambiae on the effectiveness of 
LLINs. Instead, we had to define R + and R +++ areas for the 
implementation of our study. 
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